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ABSTRACT 
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\ We present a new factorization theorem for the decay rates of P-wave states of 

heavy quarkonia. Infrared logarithms that had appeared in previous perturbative 

^ calculations of P-wave decays are absorbed into a quantity that is related to the 

Q>>^ ■ amplitude for the heavy quark and antiquark to be in a relative color-octet S-wave 



state. We predict all of the light-hadronic and electromagnetic decays rates of the 



^ • Xc and he states in terms of two phenomenological parameters. 

Ph! _ 

D . The annihilation of the heavy QQ pair in quarkonium is a short-distance 

^ ' process, occurring when the Q and Q are within 0{1/Mq) of each other. If all of the 

interactions associated with the annihilation were to occur at distance scales set 
^ ■ by Mq, then one could, because of asymptotic freedom, calculate the decay rate in 

a perturbation series in ^^(Mq). Unfortunately, quark-gluon interactions invalidate 
this simple scenario, since the associated infrared (IR) and collinear singularities 
imply the presence of long-range interactions. Nevertheless, for the annihilation 
of S-wave states a relatively simple picture emerges. If one neglects the relative 
velocity v of the Q and Q, keeping the leading term in an expansion in powers of 
v/c, then the annihilation cross section for a meson m with total spin S factors into 
a short-distance piece times a long-distance piece: 

F (m(2S+iS) ~^X) ^ Gi(m) f i (QQC'+'S) X). (1) 

The quantity fi is the short distance piece. It is the (on-shell) parton-level annihi- 
lation cross section, and asymptotic freedom allows its computation in perturbative 
QCD. The perturbation series for the annihilation of an on-shell QQ in an S wave is 
well behaved: final-state IR and collinear divergences cancel according to the KLN 
theorem, and initial-state IR divergences cancel because the meson is a color singlet. 
Gi is the long-distance piece, which contains all of the nonperturbative effects. It 
is analogous to a parton distribution, d is proportional to the probability to find 
the Q and Q at the same point: 

Gi(m)«(3/27r)|i?„,5(0)|VMQ, (2) 
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where i?ms(0) is the nonrelativistic radial wave function at the origin. 

One might guess that a factorization formula similar to Eq. (|l|) would hold 
for P-wave decays as well. However, in this case the situation is more complicated. 
Owing to the angular dependence of the wave function, the leading term in the vjc 
expansion vanishes upon integration over the angular orientation of v. Consequently, 
the annihilation cross section for P waves is suppressed by v'- jf? relative to S-wave 
annihilation. The first subleading term in the vjc expansion yields, in perturbation 
theory, a contribution that is proportional to the derivative of the wave function at 
the origin. 

Now, the first subleading term in the vjc expansion measures the color cur- 
rents of the Q and Q, rather than their color charges. Thus, for the subleading 
term, infrared divergences need not cancel, even though the quarkonium has no 
net color charge. In fact, the perturbation series for QQ annihilation on shell in 
a P wave contains a logarithmic IR divergence. This is a clear signal that long 
distance effects are present and that the use of perturbative QCD is not valid. Past 
(nonrigorous) treatments of P-wave decays have invoked the confinement scale or 
the binding energy as an IR cutoff — but with no fundamental justification. 

The structure of these IR divergences has a simple physical interpretation. 
The divergences arise when the P-wave color-singlet state converts to an S-wave 
color-octet state through the emission of a soft gluon. Then the Q and Q annihilate 
from the S-wave state. The soft-gluon transition costs a factor w^/c^. However, 
S-wave annihilation is enhanced by a factor (? jv^ relative to P-wave annihilation, so 
the color-octet process is competitive with direct color-singlet P-wave annihilation. 

By taking into account the fact that the IR divergences are associated with 
the color-octet mechanism, one can write a new factorization theorem for P-wave 
decays 

r (m(2S+ip) ^ X) - i/i(m) Fi (QQ(2^^+1P) X) + H^{m) Fg (QQf^+iS) ^ X). (3) 

The first term in Eq. corresponds to the naive (color-singlet) factorization pic- 
ture. The second term is new and gives the contribution of the color-octet mech- 
anism. The f's are the parton-level cross sections for on-shell QQ annihilation, 
except that the IR divergent part of fi is extracted and put into H^. The precise 
way in which this is done is the factorization prescription. Ei is proportional to the 
derivative of the P-wave color-singlet QQ wave function at the origin: 

i/i(m) « (9/2^)|i?;,p(0)|VAf|. (4) 

is related to the amplitude to find the Q and Q in a relative color-octet S-wave 
state. Since H% contains information about the QQg Fock state, it is not simply 
expressible in terms of the nonrelativistic QQ wave function. In perturbation theory 

where is the binding energy. The logarithm signals that Eq. (|^) is not trustworthy 
as an estimate of H^. However, by comparing it with the IR divergent parts of 



previous calculations of P-wave decayi, one can extract the f's. Both Hi and Hg have 
precise definitions in terms of operator matrix elements, so they can, in principle, be 
measured in lattice simulations, or they can simply be treated as phenomenological 
parameters. 

In general. Hi and Hs depend on the total angular momentum J and the 
total spin S of the P-wave state. However, if one describes the heavy QQ system 
in terms of a low-energy effective Lagrangian, then the spin-dependent terms are 
suppressed by powers of 1/Mq, that is, by powers of v/c. Thus, to leading order in 
u/c, we can take Hi and Hs to be independent of S and J. Corrections to this give 
terms of order v^/c^ in the decay rates, where v'^/c^ « 20% for charmonium. 

The factorization formula Eq. (^) combined with the leading-order expres- 
sions for the f's gives 

„ _ 45 r(xc2 ^ LH) ~ rjxci ^ LH) _ 1 r(xci ^ LH) 

where "LH" denotes light hadrons. Using the Particle Data Group! values for the 
branching ratios and the recent E760 datai for the Xci and Xc2 total widths, we 
find that Hi = 15.3 ± 6.6 MeV and iJg = 3.2 ± 1.4 MeV . The quoted error is the 
experimental error, which includes the uncertainty in as, combined in quadrature 
with our estimate of the theoretical uncertainty, which includes u^/c^ corrections and 
higher-order perturbative QCD corrections. Given Hi and H^ and the leading-order 
expressions for the f's, we can use Eq. @ to compute the partial widths for the 
decays of the Xco and he into light hadrons, the decay of into 7 plus light hadrons, 
and the decays of Xco and Xc2 into two 7's. There is also a simple relationship 
between the radiative decay rates of the P-statesi, which is correct to leading order 
in v/c: 

r(ipi 7 'So)/i?^(ii) ~ r(3p,, ^ 7 3Si)/ii;3(3J), j = 0, 1, 2, (7) 

where E^(ll) and E^{^J) are the energies of the 7's in the singlet and triplet decays, 
respectively. Eq. (|^ is well satisfied for the Xci and Xc2- We use it to obtain 
predictions for the radiative decay widths of the Xco and he. Our predictions for the 

Xco are 

r(xco) - (5 ± 2) MeV, B(xco-LH) = (98±l)%, 

\°) 

Bixco ^ 7 + J/^J) = (2 ± 1)%, S(xcO -> 77) = (7 ± 4) X 10"^ 

The predictions for the Xco total width and branching fraction into 7 + J/ijj differ 
significantly from the accepted values of 14 ±5 MeV and (0.66 ± 0.18)%, respectively. 
More precise data on the Xco would provide useful tests of the QCD predictions. 
Our prediction for the branching fraction of the Xc2 into two 7's is 

S(xc2 ^77) = (4±2) X 10-^ (9) 

This is somewhat smaller than the old Particle Data Group! value of (11 ±6) x lO^''. 
However, a new E760 measurement! yields a branching fraction into two 7's of 
(1.7 ±0.6) X 10^**. Our predictions for the he are 

r{hc) = (1.0 ± 0.2) MeV, B{he LH) = (52 ± 11)%, 

-B(/ic^?7c + 7) = (46±11)%, B(/ie^7 + LH) = (2±1)%. ^^^^ 



The prediction for the he total width is consistent with the upper bound of 1.1 MeV 
obtained recently by the E760 collaboration.! We predict a significant rate for he 
into 7 plus light hadrons. The hard 7 recoiling against a jet plus soft hadrons could 
be a distinctive signature for this decay. A large component of the error in all 
of these predictions arises from the theoretical uncertainty. It could be reduced 
considerably by making a complete next-to-leading-order calculation of the f 's. 

One can treat quarkonium production processes using techniques that are 
very similar to those that we have described for quarkonium decay processes. Pos- 
sible applications include photoproduction, leptoproduction, and hadron-hadron 
production. The application to production of Xc states in i?-meson decay can be 
found in Ref. ^. The nonperturbative quantities Gi and Hi, which appeared in 
quarkonium decay, appear in the color-singlet production processes as well. How- 
ever, color-octet production involves a new nonperturbative quantity H^. Whereas 
Hs is analogous to a parton distribution function, is analogous to a fragmen- 
tation function. The two quantities are related by crossing, but that relationship 
is a simple one only in lowest-order perturbation theory. Consequently, one must 
determine from experiment or extract it from a lattice calculation. 
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